A new strategy to find bound states in 
aniiarmonic oscillators 

Francisco J. Gomez and Javier Sesma 

Departamento de Fisica Teorica, 
Facultad de Ciencias, 
50009 Zaragoza, Spain. 

February 1, 2008 



Abstract 

A very simple procedure to calculate eigenenergies of quantum an- 
harmonic oscillators is presented. The method, exact but for numerical 
computations, consists merely in requiring the vanishing of the Wronskian 
of two solutions which are regular respectively at the origin and at infinity. 
The first one can be represented by its series expansion; for the second 
one, an asymptotic expansion is available. The procedure is illustrated by 
application to quartic and sextic oscillators. 

Anharmonic oscillators have been frequently used to modelize different prob- 
lems in all branches of Physics. In fact, anharmonic terms are unavoidable to 
represent a realistic potential in the vicinity of a local minimum. Besides this, 
since the seminal paper by Bender and Wu on the quartic oscillator plenty 
of mathematical questions having anharmonic oscillators as their subject have 
risen. Consequently, a great diversity of (perturbative, iterative, variational, 
etc.) methods have been proposed to find the eigenvalues of the Hamiltonian. 
Here we consider an isotropic D-dimensional anharmonic oscillator and present 
a very simple method that provides an exact quantization condition: the eigen- 
values turn out to be the zeros of an easily obtainable function of the energy. Of 
course, the determination of those zeros cannot be done without having recourse 
to numerical computation, but the accuracy of the results can be controlled. 

Let us consider the Schrodinger equation (in units fi^ /2m — 1) 

- u"{r) +V{r)u{r) = Eu{r), (1) 
with an effective (including centrifugal terms) anharmonic potential, 

2N 

V{r) = ' ^ ^ 2, A2N > 0, (2) 
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all parameters being assumed real. Such a potential may correspond either 
to a one-dimensional oscillator, with — oo < r < +00 and \u\^dr = 1, 

or to an isotropic Z3-dimensional {D > 1) oscillator, with < ?' < +00 and 
Jf^°° \u\^dr — 1. In the second case, the reduced radial wave function must 
satisfy the condition u{0) = 0. Two independent power series solutions of Eq. 
, of the form 

00 

^/(r) = ^a„r"+^ a^^O, (3) 



n=0 



with 



= ^ (1 ± ymm) , (4) 

can be immediately obtained by direct substitution. But it may occur, depend- 
ing on the values of D and that only one of such mathematical solutions, 
Urog('')j satisfies the requirement of regularity at r = 0. The normalizability of 
the wavefunction imposes a regularity at infinity that takes place only if the en- 
ergy E adopts one of the eigenvalues Ek, fc = 0, 1, 2, . . .. Asymptotic expansions 
of the wavefunction for large r, like 



u{r) ~ exp 



N+l 
P=l 



^/i„r-™, ho^O, (5) 



m=0 



can also be found by substitution. In fact, two independent solutions, u^^^{r) 
and u^'^\r), exist having asymptotic expansions of the form 0. In what follows, 
we will assume that u^^\r) represents the solution well behaved for r —>■ 00 and 
u(^'(r) a divergent one. Obviously, be written as a linear combination 

of u'^-'^' and u^^^ with coefficients T^^^ and T^^\ usually called connection factors, 
which depend on the parameters Aj of the potential and on the energy E. For 
a given set of potential parameters, the eigenvalues Ek turn out to be the zeros 
of T^^) as a function of E. 

Explicit expressions of T^^-* in terms of E can be obtained only in some 
particular cases. Nevertheless, T*^^^ can be written as the quotient of two (in- 
dependent of r) Wronskians, 

t(2) = WKog,ii''']/W[u^'\u^'^]. (6) 

As far as we are not interested in the precise value of T'^' but in their zeros as 
a function of E, the determination of the eigenvalues reduces to find the values 
of E for which the Wronskian 

WKcg,w^'^] =0. (7) 

As it can be seen, this quantization condition is most simple. Its implemen- 
tation, however, requires certain tricks, because of the formal nature of the 
expansion ((Sj. We are going to illustrate the procedure in two particular cases 
that have received considerable attention: the quartic anharmonic oscillator and 
the sextic one. 
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The quartic oscillator, corresponding to take 

T/(r) = A4r4 + A2r2+^_2r•-^ A4 > 0, (8) 

in |j2Jl, is the best known anharnionic oscillator. Its one-dimensional version 
has served as test bed for different approximate methods of solution of the 
Schrodinger equation. For the asymptotic expansions of the form © we have 
the two sets of exponents 



af^ = VA;, af^ = 0, ^ = ^2/2^^, ^(2) ^ ^ > 

and the recurrence 

2a^mh„i = {E+al)hm-i-2ai{m~l)hm-2 + iim-l){m-2)-A^2)hm-3- (10) 

The power series for Urcg and the asymptotic expansion for u^^' would allow one 
to obtain a formal expression for the Wronskian in For technical reasons, 
we prefer introducing two auxiliary functions 

v,,,{r) = exp(-4^VV3) u,,g{r), v^^\r) = exp(-aWrV3) n^'Kr), (11) 
whose Wronskian turns out to be 

W[frcg,t^''^] -CXp(-24'VV3)WKeg,u(l)]. (12) 

By substituting v^^^ by its asymptotic expansion and denoting 

wir) = exp(a^^V) Ureg(?'), (13) 
we obtain for the left hand side of H12|l 

/ oc \ 



— m— 1 

rn—Q 



+ w{r)[Y^h[l)r-^-A . (14) 



\m=0 / 

The power series expansion 

oo 

w{r) =Y,bnr'^^'' . fooT^O, (15) 

can be easily obtained as solution of the differential equation 

+ 2(4'V2 - a^p)w' + {-2A2r^ + 24'V + £;+ (aJ'V - ^-2^"')^ = 0. (16) 
The recurrence obeyed by the coefficients is 

n{n-l + 2v)b^ = 2a^^\n - 1 + v)hn-i -~ {E + {a[^^f)bn-^2 

- 24^ V - 2 + iy)bn-3 + 2A2fe„_4. (17) 
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By choosing, for instance, the starting values 

/ij,^^ = 1 and 60 = 1, (18) 

it is immediate to obtain, by means of the recurrences H1(J|) and H17|l , an explicit 
form for the asymptotic expansion (|14|l . that, with an obvious notation, would 
read 

00 

WKeg,i;(^)]^ lur'^-^+r (19) 

k—~oo 

According to p2ll , the last expression must be a formal expansion of 

exp(-24'VV3) WK-cg, (20) 

the Wronskian of Urcg and u^^^ being independent of r. The crucial point of our 
procedure is to recall the Heaviside's exponential series 



n— — 00 



introduced by Heaviside in the second volume of his Electromagnetic theory 
(London, 1899), probed by Barnes [2] to be an asymptotic expansion for ar- 
bitrary 5 and |arg(t)| < tt, and used by Naundorf in his treatment of the 
connection problem, which has provided the inspiration for our procedure. An 
asymptotic expansion for the expression (|2()|l can be immediately written. It 
becomes 

exp(-24'VV3) WKcg, u^'^] ~ Pi£i + /32f 2 + /33f3, (22) 



where 



((-2/3)aW.3)»+./3 



n— — oo 



- ((-2/3)4^V3)»+(.+l)/3 

' ^ „f-^ r(n+l + (. + l)/3) ' ^'^^ 
((-2/3)aWr3)»+(.+2)/3 
^ r(n+l + (^ + 2)/3) ' ^^^^ 

and the constants /?2 and /^s are arbitrary, with the only restriction 

WKg,u(i)] =/3i+/32 + /?3. (26) 

Identification of the right hand sides of H19|l and (|22|) allows one to obtain the 
values of (32 and /Js, that, substituted in (jSEJ, give 

WKcg, w^'^] ^77^ ( r(n + 1 + v/i) 73„+i 

((-2/3)4^V+''/3 V 

, r(n + l + (^ + l)/3) , r(n + l + (t. + 2)/3) A 

"1 rn ', 73n+2 H TT\ ] 73n+3 , {'^ ' ) 

((-2/3)4^^/3 ((-2/3)4^^)2/3 j 
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En 


E^ 
^1 


^2 







1.06036209 


3.79967303 


7.45569794 


11.64474551 


-1 


0.65765301 


2.83453620 


6.16390126 


10.03864612 


-2 


0.13778585 


1.71302790 


4.78242971 


8.33286819 


-3 


-0.59349330 


0.37766207 


3.34533567 


6.52498666 


-4 


-1.71035045 


-1.24792249 


1.94143719 


4.61294345 


-5 


-3.41014276 


-3.25067536 


0.63891956 


2.58121627 


-6 


-5.74819052 


-5.70679252 


-0.72394168 


0.37528499 


-7 


-8.67110521 


-8.66245222 


-2.54370521 


-2.11199938 


-8 


-12.13633072 


-12.13481435 


-5.12655020 


-5.01091331 


-9 


-16.12618646 


-16.12595855 


-8.44212291 


-8.41871412 


-10 


-20.63357670 


-20.63354688 


-12.37954379 


-12.37567372 



Table 1: The four lowest eigenenergies of the double- well potential H29|) . for 
several values of A2. The energies Eq and E2 correspond to even states {ly = 
in Q); El and E3 to odd ones {v = 1). 



where the index n can be chosen at will. For positive values of the subindex k, 
the coefficients of the expansion H19|l are given by 

00 

7fc = ^ /i^^ (2a^'^ bk+m -{2m + k + 2 + iy) bk+m+i) ■ (28) 

In order to compare with reliable published results, we have applied our 
procedure to the one-dimensional potential 

V{r)^r^ + A2r\ (29) 

discused by Balsa et al. Their results for A2 <Q are quoted in all posterior 
papers dealing with different approximate methods to solve the double-well 
anharmonic oscillator. Table 1 shows the four lowest eigenenergies, for several 
values of A2, obtained by requiring the cancellation of the right hand side of 
1)2 7(1 . In the computation we have used a FORTRAN program with double 
precision. Our results are in perfect coincidence with those for Eq and Ei given 
in Ref. 

The second example of the application of our strategy refers to the sextic 
anharmonic potential 

V{r)^AQr^ + Air'^ + A2r^+A^2r-^: A > 0, (30) 

which, besides having served as a model for a great variety of systems, like 
interacting anyons, molecules of ammonia, and hydrogen-bonded solids, is a well 
known example of quasi-exactly solvable (QES) potential for certain values of 
the parameters. For the asymptotic expansions (jSJ, the non- vanishing exponents 
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are now 



2 



2VA ^ SAeVA 

,(2) _ /4- ^(2) _ Aj 12) _ 3 , 4^16^2 - ^4 
V. _ i/yig Qfn — 7^=, — —77 H 7= — 



(31) 



and the recurrence 

2a4mhm = {E + a2{—2m + 5 + 2/i))/i,„_2 

+ ((m-4-/L()(m-3-/n)-^_2)/i„x-4. (32) 

As auxiliary functions we choose now 

z;reg(r) = exp(-4^VV4) u^.^ir), v^'\r) = exp(-4^VV4) u^^Hr), (33) 

whose Wronskian verifies 

W[vreg, w^'^] = exp(-af VV2) WKg, (34) 

By substituting u^^^ by its asymptotic expansion and denoting 

w{r) ^ exp(a^^VV2) Wrcg(r), (35) 
we obtain for the left hand side of H34|l 

/ 00 \ 
WKeg,^;^^)] ^ (2aW^«;(r) - w' {r)) ^ /i^^) 



r 

m=0 



^(0 E ■ (36) 



\m=0 



Since w{r) obeys the differential equation 

w" + 2{ai^\^ - J^^\)w' 

+ (-2Air^ + (3ai'^ + {a^^^ f - A^V +E- 4'^ _ A^2r-^)w = 0, (37) 

the coefficients of its series expansion 



00 



;(r) =E^"^"^'' ^0^0, (38) 



can be easily obtained by means of the recurrence 



n{n-l + 2v)hn = {-E + J^'^ (2n - 3 + 2v))bn- 



2 



+ {A2 - (4'')' - 4'^ (2" - 5 + 2v))h^.i + 2A46„_6. (39) 
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The expansion (|36|) adopts the form 

oo 

W[v..,,v^'^^ J2 72.r2'=+iW^'^ (40) 



k— — oo 

where we have denoted 



72fc = 

m=0 



oo 

^2m (Sas'^ b2k+2m - {2k + im + 2 + - fi^^'>) b2k+2m+2) ■ (41) 



The right hand side of H34|) can now be written as 

exp(-af VV2) WKcg, u^'^] - /3if i + /32£2, (42) 

with 

- ((-l/2)aWr4)»+(i+--+^''-^)/4 

' ^ „f"oo r(n + l + (l + i. + /i(i))/4) ' ^ ""^ 

- ((-l/2)4^^)"+(3W-')/4 

' ^ „f:-^ r(n + l + (3 + . + /.«)/4) ' ^^^^ 

and the constants /3i and (32 being such that 

WKcg,^^^'^] -/?i + /32. (45) 

Comparison of the right hand sides of (|40() and H42|l gives the values of Pi and 
P2 and, finaUy, 

WKg, u('^] = ^ (1) ^ n),w. ( + 1 + (1 + + /i<'^)/4) 74n 

, r(n+l + (3 + z/ + MW)/4) \ 

+ — T^^W^^ — 'H'^^ 

where the index n is any integer. 

We have looked for the values of E which make the right hand side of H4t)|l 
to become zero in the case of a potential of the form H30() , namely 

V{r) =r^ - (4s + 4 J - 2) + i(4s - l)(4s - 3) r"^. (47) 

Such a potential was first discussed by Turbiner |S] and has been also considered 
by Bender and Dunne [H] and by Finkel et al. [7] to illustrate the existence of 
a family of orthogonal polynomials associated to the quasi-exact solvability. 
Whenever the parameter J takes positive integer values, that potential becomes 
QES for any value of s. In fact, J is the number of eigenvalues and eigenfunctions 
that can be obtained exactly. In Table 2 we report the results of our procedure 
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J 


Eq 


El 


E2 


E3 


-V3/4 


5.75218468 


17.19822587 


32.16502803 


49.95708442 





4.22831744 


15.21728994 


29.75575600 


47.18056681 


0.5 


2.25118051 


12.87099068 


26.92815514 


43.93525080 


1 


0.00000000 


10.46738163 


24.05793091 


40.65160026 


1.5 


-2.56314639 


8.00423193 


21.15280932 


37.33449367 


2 


-5.46410162 


5.46410162 


18.22120991 


33.98950899 


2.5 


-8.71290067 


2.81198082 


15.27080313 


30.62297181 


3 


-12.30551201 


0.00000000 


12.30551201 


27.24194193 


3.5 


-16.22899813 


-3.02264936 


9.32087564 


23.85400371 


4 


-20.46665929 


-6.29920011 


6.29920011 


20.46665929 



Table 2: The four lowest eigenenergies of the potential with s = (2 + \/3) /4 
and different values of J. Positive integer values of J give QES potentials. The 
value J — — -\/3/4 corresponds to an effective potential + (l/2)r~^. 



by calculating the lowest eigenenergies for a particular value of s and different 
values of J. For J = 1, 2, 3, 4, the J lowest values of the energy coincide exactly 
with the zeros of the corresponding critical polynomials of Bender and Dunne. In 
these cases, the coefficients h„i vanish for m larger than a certain M , depending 
on J, and, consequently, the expansion ([SJ becomes exact and gives the quasi- 
exact eigenfunction. 

The above described method of solving the Schrodinger equation with anhar- 
monic potentials is exact, in principle, but its usefulness lies on two assumptions: 
the convergence of series like those in (|28|l and H41I) and the numerical stabil- 
ity of recurrences like l(Tn|) . lfT7|l . and in the examples considered. 
Investigations tending to justify both assumptions arc in progress. 

It is a great pleasure to dedicate this contribution to Prof. Alberto Galindo 
on occasion of his seventieth birthday. Financial support from Comision Inter- 
ministerial de Ciencia y Tecnologi'a is acknowledged. 
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